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Abstract 

We give further support for our conjecture relating eigenvalue distributions of the Kapustin- 
Willett-Yaakov matrix model in the large N limit to numbers of operators in the chiral ring 
of the corresponding supersymmetric three-dimensional gauge theory. We show that the 
relation holds for non-critical R-charges and for examples with M = 2 instead of M = 3 
supersymmetry where the bifundamental matter fields are nonchiral. We prove that, for 
non-critical R-charges, the conjecture is equivalent to a relation between the free energy of 
the gauge theory on a three sphere and the volume of a Sasaki manifold that is part of the 
moduli space of the gauge theory. We also investigate the consequences of our conjecture for 
chiral theories where the matrix model is not well understood. 
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1 Introduction 



For those interested in superconformal gauge theories in three dimensions, the matrix model 
of Kapustin, Willett, and Yaakov fl] provides a powerful tool. Using this matrix model, 
one can compute the partition function and the expectation values of supersymmetric Wil- 
son loops on a three sphere, even when the gauge theory is strongly interacting. The matrix 
model was derived through a localization procedure [2] that obscures the connection between 
matrix model quantities and microscopic degrees of freedom in the gauge theory. Given the 
success of the matrix model in post-dieting the N 3 ^ 2 large iV scaling of the free energjj^] 
of maximally supersymmetric SU(N) Yang-Mills theory at its infrared fixed point [3], it 
is a worthwhile goal to try to figure out the map between matrix model and gauge theory 
quantities in greater detail. In ref. El, we made some progress in understanding the relation 
between the eigenvalue distributions in the matrix model and the chiral ring of the super- 
symmetric gauge theory for the so-called necklace quivers, and we conjectured this relation 
would hold more generally. In this paper, we work out further examples of field theories 
that obey this conjecture. We restrict ourselves to field theories with M-theory duals of the 
Freund-Rubin type AdS± x Y, where Y is a Sasaki-Einstein space. A calculation in 11-d 
supergravity relates the field theory free energy to the volume of the internal space Y on the 
gravity side through the formula [5] 



in the matrix model and the chiral ring for the necklace quiver gauge theories. These field 
theories have M = 3 supersymmetry (SUSY), gauge group U(N) d , and associated Chern- 
Simons levels k a , a = 1, . . . , d, such that k a = 0. The matter sector consists of the 
bifundamental fields X a>a+ i and X a+l a that connect the gauge groups together into a circle 
(see figure [T]). The localization procedure [I] reduces the partition function to an integral 
over d constant N x N matrices a a , where a a is the real scalar that belongs to the same 
M = 2 multiplet as the gauge connection. In the large N limit, the matrix integral can 
be evaluated in the saddle point approximation. As was shown in [H], at the saddle point, 
the real parts of the eigenvalues of a a grow as N 1 / 2 while their imaginary parts stay of 

^y "free energy" we mean minus the logarithm of the path integral on S 3 , with an appropriate subtraction 
of UV divergences. 





where the normalization of the metric on Y used to compute Vol(Y) is R mn = 6g mn . 

Let us begin by describing the relation noticed in El between the eigenvalue distribution 
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order one as iV is taken to infinity. In addition, to leading order in N the real parts of the 
eigenvalues are the same for each gauge group. Therefore, in order to find the saddle point 
one can consider the large N expansion 



Af .V 2 Xj • n,,,, • ... . (1.2) 

As one takes iV — > oo, the Xj and y a j become dense, and one can pass to a continuum 
description by defining the distributions 

1 N 1 N 

p(x) = lim - Xj) , p(x)y a (x) = lim — V y tt jS{x - Xj) . (1.3) 

iV— >oo iv ' iv— »oo iv ' 
3=1 3=1 

The saddle point is then found by extremizing a free energy functional F[p,y a ] under the 
assumption that p is a density, namely that p{x) > and J dx p{x) = 1. It is convenient to 
enforce the latter constraint with a Lagrange multiplier p that will appear in the formulae 
presented below. In general, F[p,y a ] may be a non-local functional because the eigenvalues 
could interact with one another through long-range forces, and if this is the case the saddle 
point equations are usually hard to solve. The key insight in solving the saddle point equa- 
tions in [5] was that, luckily, in the continuum limit (1.3) the ansatz (1.2) leads to a local 



expression for F[p, y a ] due to the cancellation of long-range forces. By solving the saddle 
point equations, it was shown in [5] that the distributions p{x) and p(x)[y a (x) — yt(x)] can 
be identified for any a and b with piecewise linear functions with compact support. While 
the free energy F can be calculated by evaluating the functional F[p, y a ] on the saddle point 
configuration, it is also possible to calculate F by noticing that F[p,y a ] satisfies a virial 
theorem that gives F = 4irpN 3/2 /3 (i). 

The chiral ring of the necklace quiver gauge theories consists of gauge invariant products 
of the X a ^ a+ i and X a+ i yCL fields and monopole operators modulo superpotential and monopole 
relations. While one can define monopole operators that turn on any number of flux units 
through each U(N) gauge group, at large iV the only relevant ones are the so-called "diagonal 
monopole operators" that turn on the same number of units of flux through the diagonal 
£7(1) subgroup of each U(N) gauge group. Operators in the chiral ring therefore have an 
associated R-charge r and a (diagonal) monopole charge m. We can also introduce the 
the function ipx ab { r , m ) that counts in the same way operators that don't vanish when the 
bifundamental field X a b is set to zeron In El we found the following relation between the 



2 In our conventions, X a b transforms under the (N a ,Nb) representation of U{N a ) x U(Nb). 
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saddle point eigenvalue distribution and the chiral ring: 

d 3 tp 



dr 2 dm 
d 2 ^x ab 



drdm 



= -p{x) , (1.4a) 

m=rx//i ft 

= -p(x)[y b (x) - y a (x) + R(X ab )} . (1.4b) 

m=rx//i ft 



In other words, the matrix model eigenvalue density p(x) and the quantity p(x)[y b (x) — 
Ua(x) + R{X a b)]i which as mentioned above are linear functions of x, should be interpreted 
as derivatives of numbers of operators whose monopole charge to R-charge ratio is given by 
x/p. 

One of the goals of the current paper is to provide further evidence for the conjectures 



(1.4) in superconformal theories with gravity duals that preserve only M = 2 supersymmetry 
as opposed to the J\f = 3 SUSY of the necklace quivers studied in El. In an J\f = 2 theory, 
the U(l) R-symmetry can mix with other Abelian flavor symmetries, so the matter fields can 
have R-charges different from the canonical free-field value 1/2. The generalization of the 
Kapustin-Willett-Yaakov matrix model to non-canonical R-charges was worked out in [6j[7] . 
Furthermore, since the U(1)r symmetry can now mix with other Abelian flavor symmetries, 
it was conjectured in [6] that the correct R-symmetry in the IR can be found by extremizing 
the free energy F as a function of all trial R-charges that are consistent with the marginality 
of the superpotential. It has been seen in many examples [8 16 that this extremum is a 
maximum and that F is positive]^] 



We find that eqs. (1.4) are satisfied for more general quiver gauge theories where the 
bifundamental matter multiplets are non-chiral, meaning that they come in pairs of conjugate 
representations of the gauge group. In the first half of section [4j we examine the necklace 
quiver gauge theories, this time with an arbitrary R-charge assignment consistent with the 
marginality of the superpotential. In the second half of section [4] and appendix [Bj we 
examine theories where we add flavor (meaning M = 2 matter multiplets that transform 
in the fundamental or anti-fundamental representation of one of the gauge groups) to the 
maximally SUSY J\f = 8 theory and to the J\f = 6 ABJM theory of fl~8l. Lastly, in appendix 



B.2, we consider a theory that shares the same quiver with its (3 + l)-dimensional cousin 
that has a C 3 /Z 2 x Z 2 moduli space (see figure [3]). In all of these examples, eqs. (1.4) are 
satisfied for any choice of trial R-charges. 



Another goal of this paper is to relate the conjecture (1.4) to the observation made 
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It was suggested in [9] that F might be a good measure of the number of degrees of freedom even in 



non-supersymmetric field theories. See also 17 
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in (9 ,11 that, as checked in a number of examples, the relation (1.1 ) between the free energy 
and the volume of the internal space holds for any trial R-charges, and not just the ones that 
extremize F. That this relation)^] holds for any trial R-charges is surprising because only for 
the critical R-charges does there exist a known 11-d supergravity background AdS^ x Y. For 
non-critical R-charges, measured geometrically in terms of the volume of some corresponding 
five-cycles of Y, one can still identify a class of Sasakian metrics on Y and compute their 
volume. The volume Vol(Y) is a function of the Reeb vector of Y, which parameterizes the 
way the U(1)r symmetry sits within the isometry group of Y. We show in section [3] that 



(1.4a) holds for some choice of trial R-charges if and only if eq. (1.1) holds for the same 
choice of trial R-charges for the matter fields and a range of R-charges for the monopole 



operators. We also show an analogous result that relates (1.4b) to the volumes of five- 
dimensional sub-manifolds of Y. For a gauge invariant operator constructed from a closed 
loop of bifundamental fields X ab , it must be true that (y a — y&) = 0. Given (1.4), there 
is a geometric version of this sum that must also vanish. The last part of section 3] explains 
why. 

There are previously recognized difficulties, involving cancellation of long-range forces, 
in using the matrix model to study the large N limit of theories with chiral bifundamental 



fields |9J. We do not surmount these difficulties, but we investigate in section |5j what (1.4a) 
and (1.4b) predict for a theory with a moduli space that is a fibration over C 3 /^ (see 
figure [2]). We also study a field theory that was conjectured to be dual to AdS± x Q 2,2 ' 2 /Z k 
in appendix D (see figure [5]). 

The paper contains two heretofore unmentioned appendices. Appendix |A| proves that the 
critical R-charges maximize F for the necklace quivers. Appendix [C] reviews how to count 
gauge invariant operators for an Abelian gauge theory with a toric branch of its moduli 
space. 



2 Matrix models at non-critical R-charges 
2.1 Review of the large N limit 

To understand what it means to consider non-canonical (or non-critical) R-charges, let us 
introduce some of the ideas developed recently in refs. |6j[7j[9]. Building on the work of [I], 
refs. |6|[7) used localization to reduce the path integral of any M = 2 Chern-Simons matter 



l A similar relation between the anomaly coefficient a computed with a set of trial R-charges and the 



volume of a 5-d Sasakian space Y is known to hold in theories with AdS§ duals 19 20 . 
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on S 3 to a matrix integral. By a Chern-Simons-matter theory we mean a theory constructed 
from some number d of J\f — 2 vector multiplets with gauge groups G a (a = 1, . . . , d) and 
Chern-Simons kinetic terms i7rk a J tr A a A dv4 a + supersymmetric completion, as well as any 
number of M = 2 chiral superfields transforming in representations Ri of the total gauge 
group G = Ylt=i G a - As mentioned in the introduction, one difference between theories with 
M = 2 supersymmetry and theories with more supersymmetry is that the R-charges A» of 
the chiral fields at the IR superconformal fixed point are not fixed at the free field values 
Aj = 1/2, so the free energy will generically depend on these R-charges. In fact, it was 
proposed in [6] that a prescription for finding the correct R-charges in the IR is to calculate 
the free energy F as a function of all possible R-charge assignments consistent with the 
marginality of the superpotential and to extremize F over the set of all such assignments. 

Let us focus on the case where all gauge groups are U(N) and index the gauge groups by 
a = 1, . . . , d. Generalizing the techniques developed in j5], the authors of [9] used the saddle 
point approximation to evaluate the path integral on S 3 for a class of M = 2 Chern-Simons- 
matter theories at large N that satisfy the following five conditions: 

1. The CS levels sum to zero: Y^t=i ^ a = 0- 

2. Any matter field X transforms either in the N tt , or or (N , N&) representation for 
some a and b. 

3. The total number of fundamental fields equals the total number of anti-fundamental 
fields. 

4. For any bifundamental field X transforming in (N , N&), there exists another bifunda- 
mental field X transforming in the conjugate representation (N&, N ). 

5. For each gauge group a we have 

(R[x]-i)+ (R[x]-i)=-2. 



X m(N a ,TSS b ) Xin (N 5 ,Na) 



(2.1; 



This last condition is sufficient to guarantee the vanishing of the long-range forces on the 
eigenvalues in the saddle point approximation. Interestingly, this condition has appeared 
before in the context of superconformal (3 + l)-dimensional gauge theories. The condition 



(2.1) would imply that the NSVZ beta function of gauge group a vanishes |11|. For quiver 



gauge theories with a toric moduli space, bifundamental fields appear in exactly two terms 
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in the superpotential. Thus, if we sum (2.1) over a, we find the condition that 



(# of gauge groups) — (# of bifundamentals) + (# of superpotential terms) = . (2.2) 
In other words the quiver may give a triangulation of a torus where the faces of the trian- 



gulation are superpotential terms 21 



If these five conditions are satisfied, one can take the N — > oo limit as described in 



eqs. (1.2) and (1.3) in the introduction. The free energy is the extremum of the free energy 



functional 



F\p(x),y a (x)] = 2nN 3/2 j dxxp(x) i^k a y a (x) + A r 

\a=l 



+ 2nN 3/2 / dx \x\p( 



X in N a 



^ ,l-R{X\ 1 



+ 



TTiV 3 / 2 



2 2 " v ' V 2 T 

J dxp(x) 2 Yl + R{XWVab{x) + R[X] - l)(5y ab (x) + R[X] - 2) , 

(2.3) 



X in (N a ,N b ) 



where Sy a b(x) = y a (x) — yb{x). This formula was derived assuming the bifundamental fields 
satisfy < R[X] +Sy a b(x) < 2. Extra care must be taken when R[X] +Sy a b = or 2 because 
in these cases the discrete nature of the eigenvalues becomes important, and the equation of 



motion derived from varying (2.3) with respect to 5y a b(x) need not hold. 



Generically, the functional 2.3 has many flat directions. The following d of them play an 



important role in this paper because they correspond to changing the R-charges of the matter 
fields by linear combinations of the gauge charges with respect to the diagonal U{1) G U (N) a : 

y a (x): y a (x) y a (x) - 5 {a) , 

chiral superfield X in N a : R[X] R[X] + 5 {a) , 

chiral superfield X in N b : R[X] -»■ R[X] - 5 {b) , ( 2 .4) 

chiral superfield X in (N a , N 6 ): R[X] R[X] + 5 (a) - 5 {b) , 

A m : A m -> A m + k ^ {a) ■ 



See [9] for a more detailed discussion of these flat directions and their AdS / CFT interpreta- 
tion. 
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The A m appearing in (2.3) is the bare R-charge of the "diagonal" monopole operator T^ 1 ' . 



A monopole operator Ta 9 "- 1 turns on q a units of tr F a flux through a two-sphere surrounding 
the insertion point. Diagonal monopole operators T*" 1 ) turn on the same number m of 
tr F a flux units in each gauge group. At large N, only the diagonal monopole operators are 
important. 

We will usually impose the constraint J dx p{x) = 1 by introducing a Lagrange multiplier 
p and defining the functional 

F[p,y a ,p] = F[p,y a ) - 2nN 3 / 2 p (^J dx p{x) - . (2.5) 
This functional should be extremized with respect to p(x), y a (x), and p. 
2.2 Flavored theories 

In all gauge theories that we examine in this paper the fundamental and anti-fundamental 
fields q a and q a appear in the superpotential as 

aQa] j (2-6) 



where O a are polynomials in the bifundamental fields. It was conjectured in 22,23 that if 



this is case then the diagonal monopole operators satisfy the following OPE: 

|m| 



T (m) T (-m) = ^JJ Q a j 



(2.7) 



This OPE was conjectured in part because a parity anomaly argument shows that the 
monopole operators have gauge charges 



g a [T^]=mk a + ^Y,9a[Oo\ (2.8) 

a 

with respect to the diagonal U{1) C U(N) a , and R-charges 

R[T^} = mA m R \°^ ■ ( 2 -9) 



a 



Using the fact that each term in (2.6) must be gauge-invariant and have R-charge two, we 



have R[q a ] + R[q a ] +R[O a ] = 2 and g a [la\ + 9alq~a] + g a [Oa\ = for any a. One can use these 



S 



relations to eliminate the sum over the flavor fields in (2.3): 



F[p(x),y a (x)\ = 2vriV 3 / 2 f dx \x\ p(x) ( R[T^\ + ^2y a (x)g a [T^ 
+ j dxp{xf Yl + R[XWy ab {x) + R[X] - l)(5y ab (x) + R[X] - 2) . 



X in (N a ,N 6 ) 



(2.10) 



3 An Equivalent form of our conjecture 

3.1 Eigenvalue density and volumes of Sasakian spaces 



We can relate the conjecture (1.4) to the observation that eq. (1.1) holds for any trial R- 
charges. In particular, we prove the following result: In a CS-matter theory dual to AdS^ x Y 
fix a set of matter R-charges R[X] and a bare monopole charge A m so that the conformal 
dimensions of all gauge-invariant operators satisfy the unitarity bound. Let p(x), p, and 
ip(r,m) be as defined in the introduction, and let's assume p{x) has compact support. The 
following two statements are equivalent: 



A. The conjecture (1.4a) holds for the given R-charges R[X) and bare monopole charge A m . 



B. For any 5 in a small enough neighborhood of zero, we have 

lim fULJL. = voi(y, 6) , (3.1) 

n^oo 27F 2 K ' 

where the free energy F of the CS-matter theory and the volume Vol(Y, S) of the internal 
space Y are both computed assuming that the matter R-charges are R[X] and the bare 
monopole charge is A m + 5. 



For notational convenience, let's denote the LHS of eq. (3.1) by Vol m (Y, 6) and let's 
introduce the rescaled matrix model quantities: 

x = - , p(x) = P< ^- , y a (x) = y a (x) . (3.2) 
p p 
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The equivalence between (|A]) and (|B]) follows from the following two equations: 

4 fdx-^ . 
24 7 {l + x5Y 



Vol m (Y,6) = - I dx - Hy > , (3.3) 



\r uv x\ ^ I a- lim r ^oo^ (2,1) (r, rx)/r 



which we prove in sections 3.2 and 3.3, respectively. 

Assuming the eqs. (3.3) and (3.4) to be true, it is clear that the statement (|A]) implies 
([Bp. That (|B]) implies (|A]) follows from the fact that knowing Vol m (Y, 8) for 8 in a small 
neighborhood of zero, one can reconstruct p(x), and analogously, from Vo\(Y,8) one can 
reconstruct lim^oo ^>( 2,1 )(r, rx)/r. Indeed, one can extend Vol m (Y, 8) to any complex 8 as an 
analytic function with singularities. We assume that p is supported on x + ] for some X- < 



< x + . We see from eq. (3.3) that the integral converges absolutely if 5 G (—l/x + , —1/xJ) 
or 5 ^ K, so Vol m (Y, 5) can only have singularities on (— oo, — l/x + ] U [— oo). 

To relate the singularities of Vol m (Y, 5) to p(x) we can perform two integrations by parts 



in (3.12) 



for any 5 G C\ ((— oo, — l/x + ] U [— l/£_,oo)). Generically, eq. (3.5) shows that Vol m (Y, 5) 
has two branch cuts, one on (— oo, —l/x + ] and one on [— l/s_, oo). From the discontinuities 
of Vol m (Y, 5) one can read off p /; (— 1/5). Simple poles of Vol m (Y, S) at 5 = — 1/x' correspond 
to contributions to p"(x) proportional to 5(x — x'); second order poles of Vol m (Y, 8) at 
5 = —1/x' correspond to 5'(x — x'), etc. From the singularities of the analytic continuation 
of Vol m (Y, S) one can therefore reconstruct uniquely p"(x), and hence p(x), and similarly for 
Vol(Y, S) and \im T ^ OQ 'tp( 2,lS> (r,rx)/r. If Vol m (Y, 8) and Vol(Y, 8) agree on an open set, then 
Q holds. 

In our examples, Vol m (Y, 8) is a rational function of 8 with poles of order at most three, 
so p(x) is piecewise linear and it may have delta-functions. From the location and residues 
of the poles one can first reconstruct p"(x), and then p{x) by integrating p"(x) twice. To 
perform this reconstruction starting with Vo\ m (Y,8), one first decomposes Vo\ m (Y,8) into 
partial fractions, and then identifies the terms in p"(x) that give those partial fractions: if, 
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for example, 



Vol m (F,5) 



7T 



7T 



E 



(3.6) 



for some Xi, then 



p"{x) 



(3.7) 



3.2 Matrix model dependence on 5 



In this subsection we prove the result (3.3). As we have seen in the previous section, the 
matrix model generally takes the form 



(3i 



F[p,y a ,lA = J dx p{x) 2 f{y a {x)) - J dx p(x)V(x,y a (x)) 

+ M*./.| I | pH ^1-M^(/ if( „- 1 ,, 



for some functions / and V. While the explicit form of these function is given in (2.10) 



their precise form doesn't matter. The only property of V that we will use is that it is 
homogeneous of degree one in x, namely V(Xx,y a (x)) = XV(x,y a (x)) for any A > 0. With 



the rescaling (3.2), one can write F as 



F[p, y a , p] = -2nN 3/2 p + p 3 J dx 



x 



-^f{ya{x)) 



X 



X 



X 

1 

l.rl 



(3.9) 



The rescaling (3.2) is useful because now the equations of motion for p and y a do not 



involve p. One can first solve these equations, and then p can be found by integrating p: 
the normalization condition J dx p(x) = 1 becomes 



dx p(x) 



(3.10) 



We now see that if we extremized (3.9) in the case where the monople R-charges were 



-Rp^ 1 )], we could obtain the saddle point when they are RlT^ 1 ^] ±5^ ±lS> through the trans- 



it 



formation: 

ps(x s ) p{x) 1 1 



_|_ ^(sgni) 



x 5 x x 5 x (3.11' 



Indeed, the equations of motion for p and y a are obtained by extremizing the expression in 



the square brackets in (3.9), and this expression is invariant under (3.11). Given that p has 



compact support, the transformations (3.11) make sense only when 6^' are small enough. 



For simplicity, from now on let's restrict ourselves to the case even 
though one can make similar arguments for the case where 5^ +1 ' ) and are arbitrary or 
satisfy a different relation. In 4 , we showed that F = AttN 3 ^ 2 p/3, which implies that 



VoUY.S) = ^ = di si>s {i s ) = £J dij^^. (3.12) 
3.3 Operator counting dependence on 5 



We now prove the result (3.4). Let A be the chiral ring associated to the superconformal 
field theory dual to AdS^ x Y in the Abelian case N = 1. A is also a vector space over C 
that is graded by the R-charge and monopole charge, meaning that one can define a basis of 
operators with well-defined R-charge and monopole charge. Let A m%r be the vector subspace 
of elements of A with monopole charge m and R-charge r. We introduce the Hilbert-Poincare 
series 

f(t, u) = J2 dim(A m , r )fu m . ^ 13 ) 

m,r 

Since the Abelian moduli space of the gauge theory is the Calabi-Yau cone over Y one can 
view the operators in the chiral ring as holomorphic functions on this cone. Martelli, Sparks, 
and Yau [24] show that 

Vol(F,5) = ^lim(l-t) 4 /(t,t 5 ). (3-14) 
48 t->i 

One can compute the Hilbert-Poincare series for Y in terms of ip(r,m), the number of 
operators with R-charge at most r and monopole charge at most m. Approximating ip by a 
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continuous function of homogeneous degree four, the definition (3.13) gives 



/(*,«)« drdmij ( - 1 ' 1 \r,m)t r u m . 



Since 1 — tm — hit for t ~ 1, we can use (3.14) and (3.15) to write Vol(Y,5) as 



Vol(y,5) « ^(lnt) 4 / drdm^ 1 ' 1 \r,m)t r+mS 
48 



7T 

'48' 



= - — (hit) 3 / drdxri> (2 ' l \r,rx)f {l+&5) 
48 J 

71-4 /" a- V' ( r > rx)/r 

= 24 y (i+H) 3 ' 

where in the second line we integrated by parts once, and in the third line we defined m = rx. 



(hit) 3 / drdm^ l \r,m)f +m5 



(3.15) 



(3.16) 



3.4 Matrix model and volumes of five-cycles 

For any gauge invariant operator X, we should have [25] R[X] = ix Vol(Ex)/6 Vol(F), where 
by T, x we denoted the 5-d submanifold of Y defined by the equation X = 0. Using Vol(F) = 
(7r 4 /24) J dxp(x), one can rewrite this equation as 



Vol(S 



7T 



X, 



dx p(x)R[X] . 



(3.17) 



For an operator X that is not gauge invariant, such as a bifundamental field that transforms 



in (N a ,N&), R[X] is not invariant under baryonic symmetries (2.4), but R[X] + y a (x) — yb(x) 
is. So we suspect that 



Vol(E 



7T" 



X, 



dxp(x)(R[X] + y a (x) - y b {x)) 



(3.18) 



We can think of this relation as a conjecture and prove the following result: If X is a chiral 
operator transforming in (N Q , N&), then for 5 in a neighborhood of zero, let 



7T" 



Vol m (Ex,5) = y / rfx,5p ( 5(x ( 5)( J R[X] +y a ,«5(£<5) ~ VbA^s)) ■ 



(3.19) 



The following two statements are equivalent: 
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I. The conjecture (1.4b) holds for the given R-charges R[X] and bare monopole charge 



II. For any S in a small enough neighborhood of zero, we have 



Vol m (E x ,5)=Vol(E x ,<5) 



(3.20) 



where the volume Vol (Ex, 5) is computed with the induced Sasakian metric on Y that 
corresponds to the matter R-charges R[X] and the bare monopole charge A m + 5. 

The proof of this result is similar to that of the equivalence between (|A| and (IBl we 



discussed above, so we skip most of the details. Using (3.11), one can check that 



Vol m (S x ,5) 



7T" 



dx 



p(x)(R[X] + y a (x) - y b (x)) 



(3.21; 



4 J (1 +x5) 3 

Defining fx(t,u) to be the Hilbert-Poincare series for the ring of chiral operators obtained 



from the chiral ring by setting X = 0, and using the Martelli, Sparks, and Yau result 24 

r3 



Vol(£x,5) = 
one can show as in section 13.31 that 

Vol(V x ,6) = 4 



7T" 



lim(l-t) d /x(M d 



(3.22) 



dx- 



lim r 



[1 + x5) 3 



(3.23) 



Here, ipx{f, m) denotes the number of chiral ring operators with X = 0, R-charge at most 
r, and monopole charge at most m, and can be approximated by a smooth function of 



homogeneous degree three. By an argument analogous to the one in section 3A it follows 
that the statements (IT1) and (JTT|) are equivalent. 



3.5 A Consistency condition 

Note that gauge invariant operators in the quiver that have no monopole charge are con- 
structed from closed paths of bifundamental fields O a . A consequence of our conjecture 
( 1.4b ) is then that for a gauge invariant operator X = Y[ a O a with no monopole charge the 
following sum vanishes: 



E 



drdm 



-R[O 



dr 2 dm 







(3.24) 
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Of course J2 a R[@a] = -R[-^]> an d we can simplify this expression: 

d 3 tfj 



E 



Or dm 



R[X] 



dr 2 dm 



(3.25) 



We would like to show why (3.25) must hold from geometric considerations aloner 



The number of gauge invariant operators of fixed R-charge r and monopole charge m 
that do not contain the operator X is approximately 



(r, m) w (r, m) - (r - R[x] , m) « RlX}^ 2 ^ (r, m) 



(3.26) 



when r 3> is large. 



We can also use (3.22 ) to express the operator counts in terms of volumes. For coordinates 



x and y on a compact space, the volume of the set of points xy = is the union of the set 
of points where x = with the set of points where y = 0. That the volumes are additive 
implies the identity 

Vo1(£ x ,5) = ]TVo1(£ 0q ,5) . (3.27) 



From this identity at large r and the results in the earlier part of this section, we have 



x (r, m 



) = £* 



(3.28) 



Combining (3.26) with (3.28) yields (3.25). 



4 Theories with non-chiral bifundamental fields 

4.1 J\f = 2 deformations of the necklace quivers and matrix model 

Our first field theory example consists of deformations of the necklace quiver gauge theories 
whose (undeformed) matrix models were also studied in pip). In M = 2 notation, the field 
content of the necklace quiver theories consists of d vector multiplets with Chern-Simons 
kinetic terms and coefficients k a , and chiral multiplets A a and B a connecting the gauge 

5 For a Calabi-Yau four fold, at large r and m we approximate ip(r,m) and tpo a {r,m) by homogeneous 
polynomials of degree four and three respectively. 
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Figure 1: A necklace quiver gauge theory where the gauge sector consists of d U(N) gauge 
groups with Chern-Simons coefficients k a and the matter content consists of the bifunda- 
mental fields A a and B a . 

groups into a necklace (see figure [T]). The superpotential 



d 

W = Y,TT ^(B a+1 A a+1 - A a B a ) 



(4.1) 



0=1 



preserves M = 3 supersymmetry. For any given k a satisfying Yl a =i ka = 0, the field theory 
is dual to AdSi x Y where Y is a tri-Sasakian space, which is by definition the base of a 



hyperkahler cone 26 

While N = 3 SUSY restricts the R-charges of A a and B a to be 1/2, in this section 
we examine what happens if we make more general R-charge assignments for the A a and 
B a fields that break M = 3 down to = 2. These R-charge assignments are required to 



preserve the marginality of the superpotential (4.1). This condition implies that for generic 



values of the CS levels k a , namely if there are no cancellations between the various terms in 
(4.1), we must have R[B a ] = 1 — i2L4 a ]. The matrix model free energy functional is in this 
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case 



F[p,y a ,ij] = 2vriV 3/2 / dx px QaSy a + 2vr A m iV 3/2 / dx px 

J a=l J 

- TrN 3 / 2 fdxp 2 J2 ~ R[M) {$Va + R[Ba}) - 2<kN*' 2 [Ji ( f dx P - 1 

J a=l 



(4.2) 



where 5y a = y a ~\ — y a , k a = q a +i ~ Qa- As per the discussion after eq. (2.3), the equations of 



motion for 5y a following from (4.2) hold only when —R[B a ] < 5y a < R[A a ]. It is possible to 
have Sy a = R[A a ] or 5y a = —R[B a ], but in that case we should not impose the equation of 
motion for that particular 8y a . 



Using the rescalings (3.2), one can write the solution of the equations of motion following 



from (4.2) as 



p{x) = s L (x) - s s (x) , 5y a (x) 



R[A a ] - R[B a ] 1 \s L (x) + q a x\ - \s s {x) + q a x\ 



2 2 s L (x) - s s (x) 

where Sl(x) > ss(x) are the two solutions of the equation 



(4.3) 



s(x)d + xc 2 + 22 l s (^) + = 2 • 



(4.4) 



a=l 



with 



ci = ( R l A -} - R i B *]) > c 2 = 2A m + ^ <?a (R[A a ] - R[B a }) . 



(4.5) 



a=l 



a=l 



The constraint imposed by varying F with respect to p is f dx p(x) = l/p 2 - 

We have encountered a solution of this type in |4] in the case where -R[A a ] = R[B a ] = 1/2 
and A m = 0. As in [4], one can think of eq. (4.4) as defining the boundary of a polygon 



V = <(x,s) eR 2 : sci + xc 2 + ^2\s + xq a \ < 2 



(4.6) 



0=1 



The quantity p(x) = sl(x) — ss{x) can then be interpreted as the thickness of a constant x 
slice V x through this polygon, p(x) = Length("P £ ). Consequently, J dx p{x) = Area("P) and 



Vol m (F) 



4 4 
71 71 



24p 2 24 



Area(P) 



(4.7) 
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See appendix [A] for a proof that the J\f = 3 R-charge assignments minimize Vol m (Y) or, 
equivalently, maximize F. Just like p(x), the quantities p(x)Sy a (x) can also be given geo- 
metrical interpretations: 

p{x){5y a {x ) + = Length {V & n {s + q a x > 0}) , 

(4.8) 

p(x)(-5y a (x) + R[A a ]) = Length (V £ n {s + q a x < 0}) . 
The equations above were written in a way that makes manifest the invariance under 



the flat directions exhibited in (2.4). Indeed, while in writing the free energy functional 



(4.2) we assumed i2[A a ] and A m to be independent, we see that the eigenvalue density p(x) 



and the quantities appearing on the LHS of (4.8) depend non-trivially only on the linear 



combinations c\ and c 2 that were defined in (4.5). These are the only linear combinations 



of -RL4J and A m that are invariant under all symmetries in (2.4). The reason why we were 
able to find two such linear combinations at all is that the spaces Y have generically two 
£7(1) isometries that commute with U(1)r. 

4.2 Operator counting for necklace quivers 

We now relate the matrix model quantities p(x) and p(x)5y a from the previous section to 
numbers of operators in the chiral ring of the gauge theory when N — 1. In [4] we provided 
such a relation in the case i?[A ] = R[B a ] = 1/2 and A m = 0, and the argument presented 
in that paper holds, with minor modifications, for the more general R-charge assignments 
considered in this paper. As explained in [I], gauge invariant operators can be constructed 
out of the bifundamental fields A a and B a and the diagonal monopole operators T^ m \ and 
they are 

0(m, s,i,j) = T ^cr i+s C™ q2+s ■ ■ ■ C a nqd+s (A 1 B 1 ) i (A 2 B 2 y , 

A mqa+s iimq a + s >0 (4.9) 

B -mq a -s if + s < o 



(jrnq a - 



The labels m and s run over all integers, while i and j should be nonnegative integers. 

Let tp(r,m) (i/) (r,m)) be the number of operators 0(m,s,i,j) (0(m,s, 0,0)) with R- 
charge at most r and monopole charge at most m. In [i] we showed that at large r and m we 
have ip( 2 '°\r,m) m ipo(r,m). This relation holds for the more general R-charge assignments 
too because the only assumption needed to prove it was i2L4ii?i] = i2[A 2 5 2 ] = 1, which we 
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still assume. A simple computation yields 



R[0{m, s, 0, 0)] = mA rn + £ R[C k ] \ mq a + s\ 



sci + mc 2 + \s + mq a \ 



a=l 



(4.10) 



where C\ and c 2 are as defined in (4.5). Using this formula one can check, as in [4], that 
^( 2 ' 1 )(r, rx)/r ~ ^q 0,1 '' (r, rx) /r is indeed given by the length of the slice V& through V. We 



have therefore verified explicitly eq. (1.4a) for the necklace quivers at non-critical R-charges. 



Let ipx a ( r ! m ) be the number of chiral operators with R-charge at most r and monopole 
charge at most m that are nonzero when X a = 0. As in ji , we have that ip^\r, m) equals 
the number of operators of the form 0{m, s, 0, 0) with R-charge at most r and monopole 
charge at most m with the extra constraint that mq a + s < if X a = A a and mq a + s > if 
X a = B a . As argued in 4], these extra constraints imply that when r is large i^x ( r J rx)/r 
is given by the length of the intersection between the slice V± and the half-plane s + q a x > 
if X a = B a or s + q a x < if X a = A a . Comparing with eq. (4.8) we see that the necklace 



quivers at arbitrary R-charges also obey our second conjecture (1.4b). 



4.3 Flavored necklace quivers 

The discussion in the previous two subsections can be generalized by including flavor fields 
that interact with the existing matter fields through the superpotential 



a=l 



.3=1 3=1 



(a) 



(4.11) 



Given that the A a transform in (N a _!,N a ) and the B a transform in the conjugate repre- 
sentation (N a _i,N a ), for eq. (4.11) to make sense we must take gj , gj a \ Qj, and Qj to 



transform in N a _i, N a , N a _i, and N , respectively. 

We discussed a superpotential of this form at the end of section [2j where we found that the 
effect of including the flavor fields was that the CS levels k a and A m of the unflavored model 
were replaced by (sgn£)<7 [T( sgna: -'] and (sgnx)i?[T ( - sgn:1 '- ) ], respectively. Eqs. (2.8) and (2.9) 
applied to our flavored necklace quivers give 



sgnx 



m a - n a+l +m a+1 ) , 



k a ^(sgnx)g a [T^] = k a + 
A m -> (sgnx)R[T^} = A m + ^ E K#[A a ] + m a R[B a \) 



(4.12) 
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From k a = q a+ \ — q a we further have 



sgnx 

q a -> q a 7T~ i n o 



(4.13) 



We believe that all the formulas presented in the previous two subsections continue to hold 
for the flavored theory if one makes the above three replacements. In particular, the relation 



between the matrix model quantities and operator counting we conjectured in eq. (1.4) 



continues to hold, and the volume of the 7-d space Y is still proportional to the area of a 



polygon V of the type (4.6). 



4.4 Flavored J\f = 8 theory and its matrix model 



We broaden our scope of examples and verify (1.4) for maximally supersymmetric Yang-Mills 



theory to which we add flavor. The theory has one gauge group and three adjoint fields X; t 
i = 1, 2, 3 coupled to n\ + n 2 + n 3 pairs of fundamental fields through the superpotential 



W ~ tr 



ni ri2 nz 

Xl [x 2 , x 3] + y: # } *i<f + E *f «f + E ^ x ^. 



ft) 

3 



3=1 



3=1 



3=1 



(4.14) 



The corresponding matrix model was solved in [9] in the large iV limit. We review their 
solution for p(x). In the next subsection, we will compare p(x) with the distribution of 



operators in the chiral ring and show that (1.4a) holds. In this case, eq. (2.7) takes the 



form T^T^ 1 ) = X^X^X™ 3 22,23 . To keep the notation concise, we define A; = R[Xi] 



A = RIT^ 1 )}, and A = R[T^~% The matrix model free energy functional is then 



p[p] = vriV 3 / 2 



dxp (AiA 2 A 3 p+ (A + A) \x\ + (A - A)x) - 2p [ dxp-1 



(4.15) 



As before, we define the hatted quantities (3.2). The eigenvalue density p(x) is 



p{x) 



AiA 2 a 3 if — 2; < a; < . 
otherwise , 



(4.16) 



which agrees with (4.8) of (9). 



20 



4.5 Operator counting in flavored M = 8 theory 



The gauge- invariant operators built out of diagonal monopole operators and adjoint fields 
in this theory are trlT^ X^ 1 X% 2 X£ 3 ]. The R-charges of these operators are 



R[T {m) X^X^ 2 X^} 



mA + Ym=i a i^i m>0, 



(4.17) 



Let ip(r,m) be the number of operators with R-charge smaller than r and monopole 
charge smaller than m. To match with p(x), we want to calculate d 3 ip/dr 2 dm at large 
r. It is easiest to start by calculating the derivative dip /dm which equals the number of 
operators with R-charge smaller than r and monopole charge equal to m. For m > 0, at 
large r the number of operators tr[T^ X^ 1 X^ 2 X^ 3 ] is approximately equal to the volume 
of a tetrahedron with sides of length (r — mA)/A$; similarly, for m < 0, the number of 
operators is equal to the volume of a tetrahedron with sides of length (r + mA) / A*. We thus 
have 



dip 
dm 



6A 1 A 2 A 3 11 U ^ ^ A ' 
6A1A2A3 11 A 

otherwise . 



(4.18) 



Taking two derivatives with respect to r, we find agreement with (4.16) and confirmation of 



the conjecture (1.4a). 



4.6 Other examples 

We presented flavored M = 8 in the main text because of its simplicity. One disadvantage 
of this example is that it possesses a single U (N) factor and so we could not compute 



a 5y and check (1.4d). To remedy this problem, in appendix [B] we consider two more 
complicated examples. The first of these is ABJM Chern-Simons theory (a theory with two 
gauge groups) [18) to which we add flavor. The second example has four gauge groups (see 
figure [3]). When a four-dimensional gauge theory has the field content of this second example, 
the Abelian moduli space is a Z2 x Z2 orbifold of C 3 . Thus, with some abuse of notation, 
we refer to this second example as the Z 2 x Z 2 orbifold theory. 



The verification of (1.4) requires on the one hand calculating p(x) and 5y(x) using the 



large iV limit of the matrix model (2.3) and on the other counting operators in the chiral 
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ring. We have two methods at our disposal for this counting. One may count the operators 
directly as we did above. Because the moduli space is toric for these last three examples, 
the direct approach has some generic features which we review in appendix [Cj In section 
[3j we presented an indirect counting method that involved calculating Vol(Y, 6) (3.16) and 
Vol(£x,<5) (3.23) as a function of p(x) and Sy(x). 



5 Theories with chiral bifundamental fields 



5.1 Noncancellation of long-range forces 



As noted in j9j, the functional (2.3) does not appear to describe the large N limit of gauge 



theories with chiral bifundamental fields. To derive (2.3), it was assumed that the long-range 



forces on the eigenvalues cancel. But for theories with chiral bifundamentals, there is no such 
cancellation. 

The long-range forces at issue come from the interactions between the eigenvalues, both 
within a vector multiplet and between vector multiplets connected by a bifundamental field 
X ab [91: 



^ = £c th7r(Aj a) -Af) 



L inter 



77.(6,0) 
L inter 



E 

3 

E 



3 L 



R[X ab ] - 1 - i(A? } - A^) 



coth tt ( \\ b) -Xf- z(l -R[X ab ]) 



(5.i; 



cothvrfAf-Af + i(l-R{X ba }) 



If 



A (a) _ A (6) 



^> 1, then we may approximate coth a; ~ sgnRex. The long-range forces are 



the forces (5.1) with coth replaced by sgnRe. For theories with non-chiral bifundamentals 



and equal ranks, the long-range forces cancel out when Re A^ = Re X\ u> for all i, a, b and 



(|2.1|) is satisfied. In general, the long-range forces on A^ cancel out only when 



J2( R i x *b] - 1 + Vbj) + ^(RiXb*) - 1 - y b ,j) = -2 . 



(5.2) 



Thus the free energy functional (2.3) is correct for theories with chiral bifundamentals only 



if the y a (x) satisfy some constraints. 
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5.2 Operator counting for the C 3 /Z3 theory 

To investigate what the matrix model for a chiral theory should give in the large N limit, we 
study the U(N) 3 Chern-Simons theory described by the quiver in figure [2j Let the Chern- 
Simons coefficients be (k±, k 2 , k 3 ) such that k\ + k 2 + k 3 = 0. We will assume k± > 0, k 2 < 0, 
k 3 < 0. The moduli space is a Kahler quotient of C 5 with weights + fc_), |(A;+ + 

k-), \{k + + — —k-), where we define k- = k\ — k 2 and fc + = fci — k 3 . 
There is a superpotential of the form 



W ~tr[6 ijk A 31}k A 23d A 12 ,i} , (5.3) 
and a monopole relation T^T^ = 1. We let i^Ay,].] = A*, i?[A iii2 ] = A y , R[A ijt3 ] = A z , 



with A x + A v + A 2 



2 as any other choice of R-charges may be transformed into this choice 

: -R[T(-V] = A. 



by a transformation of the form (2.4). We denote R[T^ 



The gauge invariant operators have the form 



r (m) ni](V),-: 



n i(i+l),J 



(5.4) 



To be gauge invariant, for m > we must impose Y2j n i2,j = + s ' 2~^j n 23j = + 
mk 2 + s and J2j n 3i,j = s an d for m < we must impose . n 23 j = mfc 2 + s, ^ ■ 7131J = 
m/c2 + mk 3 + s, and X}j n i2,j — s - Given the R-charge assignments, it is convenient to 
introduce rij = J2i n i(i+i),j- Each gauge invariant operator corresponds to a quadruple 
(ni, n 2 , n 3 , m) such that . rij = m/c sgn ( m ) + 3s and m is bounded between — . rij/k_ and 

Given the description of the gauge invariant operators, it is now a straightforward task 
to count them by either the direct method described in appendix [C] or the indirect method 
described in section [3] For A x > A y > A z a piecewise expression for p(x) is: 



p(x) 



0, 

l+(fe_A z — A)x 
3A z (A x -A z )(A y -A z ) 



(A, 



-A z )(l-A£)-A !/ A z fc_x 



3(A a ,-A w )(A a ,-A JS )A, ( A a 
1-Af 



(A 3: -A i/ -A 2 )(l-A£)+A i) A z fc + £ 
3(A ;c -A y )(A ;c -A z )AyA z 
l-(fc + A z +A)£ 



3A z (A a; -A z )(A !/ -A z ) 



0. 



x < 



fc_A z -A ' 



fc_A z 



k-A x -A 
1 



<X< 
<X< 
< X < 



1 


k-Ay 


-A 


1 






-A 


1 





fc+A y +A 
1 

fc + A z +A 



< T < 
fc+A^+A fc+Ay+A ' 

1 < r < 

•' - fc+A z +A ' 



fc+A^+A ' 
1 



< X . 



(5.5) 
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We note three odd things about (5.5): 1) If A x = A y = A z , p has a delta function at 



k-A x 



and 



fc+A^+A' 



2) In contrast to nonchiral examples, p(x) while still piecewise linear 



is no longer a convex function of x. 3) The matrix model (2.3) gives the same result for p in 
the central region despite the fact that the long range forces do not cancel. (In other regions 
and for 5y a b, the matrix model results are different.) 

Now we set ^° zero - The nonzero operators are those with no A^ i+ i^i fields. As a 
piecewise function: 



-y 2 (x) = < 



-A,, 

A y A z (l+(k-A x -A)x) 
(A x -Ay-A z )(l-Ax)-A y A z k-x ' 

0, 

A y A z (l-(k+A x +A)x) 
(A x -Ay-A z )(l-Ax)+A y A z k + x ' 

-A 2 , 



1 



< x < 



fc_ A z -A 

■ , } A < X < 



1 



fc-A,-A 
1 



<X < 



1 


k-A v 


-A 


1 




k-A x 


-A 


1 





fc+A^+A 



1 



fc+A„+A 



< X < 

< X < 



k + A x +A 
1 



k + A y +A 



1 



fc + A z +A • 



(5.6) 



Finally, we set ^3^1 to zero. The nonzero operators are those with no A^+i^i's, and 
those with m > 0, n x + n y + n z = k + m. As a piecewise function: 



yi{x) - y3[x) 



-A,, 

A y A z (l+(fc_A^-A)x) 
(A 2 ,-A !/ -A z )(l-A£)-A a A z fc_ : r ' 

0, 

A y A z (l-(k+A x +A)x) 
2 A, 



k-A z - 


A 


1 




k-Ay- 


A 


1 




k-A x - 


A 


1 


< 


fc+A^+A 


1 


< 



< x < 
<x< 
<x< 



k-A x -A ' 
1 



<X< 



k + A x +A ' 
1 



(5.7) 



fc + A z +A ' 



The result for i)2{x) — yi(x) follows by taking the difference of (5.6) and (5.7). We have 
checked that the operator counts where we set each of the remaining seven bifundamental 
fields to zero in turn yield the same results for the differences in the i/'s. 



5.3 Particular case: the cone over M ' ' /Z& 



Consider the case where the internal space Y is M 1 ' 1,1 /Zfc. It was proposed in 27-29 that 
the dual field theory is the one in figure [2] with CS levels k\ = 2k and k 2 = k% = —k, so 
k + = k- = 3k. As a function of the trial R-charges, the volume of Y is 



Voi(y) 



7T 
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dx p(x) 



3fc 3 7r 4 (A 2 + 9k 2 (A x A y + A X A Z + A y A z ) 
8 (9fc 2 A2 - A 2 ) (9k 2 A 2 y - A 2 ) (9k 2 A 2 z - A 5 



(5i 
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Figure 2: The quiver for the C 3 /Z 3 theory. When the CS levels are (2k, —k, —k) this field 
theory is believed to be dual to AdS^ x M 1,1 ' 1 /Zfc. 



Under the constraint A^+A^+A^ = 2, this expression is maximized for A x = A y = A z = 2/3 



and A = 0, and the maximum is 97r 4 /(128fc), which is the volume of M 1,1 ' 1 /Z^ 25 
For the critical R-charges, our predicted eigenvalue density is 



p(x) (y 3 (x) 
p(x) (y 2 (x) 



p(x) 

m(x)) 
m(x)) 



8 \2k 1 



3 . ( „ 1 
— -6 [ x + — 
16k \ 2k 



9 A _1_ 
32k [ X+ 2k 



9 
32k 



5 x 



1 

2k 



3 J, 1 
16**1* + 2* 



16k 
3 



-5 ( x 



1 

2fc 
1 



8A; V 2k 



The volumes of the five-cycles corresponding to the bifundamental fields are 

- 3 r f 2\ 3vr 3 



voi(s A23a ; 

Vol(S Al2 J 
Vol(S A31 J 



dxp(x) y 3 (x) -y 2 (x) 



— / dxp(x) (yi(x) - y 3 (x) + - 



16k 

217T 3 

~6Ak 
21tt 3 
64k 



(5.9) 



(5.10) 



Let us understand how these volumes are related to the volumes of the divisors computed 
in [25]. The cone over M 1 ' 1,1 is a Kahler quotient of C 5 by a U(l) that acts with weights 
(2,2,2,-3,-3) on the coordinates (u\, u 2 , W3, v 1, v 2 ) parameterizing C 5 . The orbifold 
used to produce the quiver in figure [2] acts by the identification (v\, v 2 ) ~ (v\e 2m ^ k , v 2 e~ 27r ^ k ) 
leaving the Ui coordinates untouched. It is natural to identify A 23a with u a , Ax 2>a with u a Vi, 
and A 3 i a with u a v 2 . Using the explicit metric on M 1,1 ' 1 the authors of [25] calculated the 
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volumes of the five-cycles corresponding to either u a = or v b = in M ' ' to be 

Vol(E u J = ^, Vol(E„J = ^. (5.11) 

We see that these equations are consistent with ( |5.10[ ): we have k Vo1(£a 23i J = Vol(E Ua ) as 
well as k Vol(E Al2i J = Vol(E u J + Vol(E B1 ) and k Vol(E A31 J = Vol(E u J + Vol(E V2 ). The 
factor of in these formulas comes from the fact that the cycles whose volumes are given in 



( pTTot belong to a Z k orbifold of M 1 ' 1 ' 1 . 

For those interested in another simple example of a theory with chiral bifundamental 
fields, we describe our predictions for a theory with the cone over Q 2 > 2 > 2 as its Abelian 



moduli space in appendix D 



5.4 Missing operators 

There is a difference between the matrix model and operator counting that manifests itself 
in chiral theories. The matrix model depends explicitly on the bifundamental fields, and a 
Sy saturates when it reaches minus the R-charge of a bifundamental field. In the absence 
of flavors, the saturation of the Sy is responsible for all of the corners in p and p5y. In the 
C 3 /Z 3 example, p has a corner at x — fc+A X + A - We might expect that there exists some 
bifundamental field so that 5y + R[Aij^] becomes zero at x = fc Aa , +A , or equivalently 
that ( r > ™) becomes zero at x — k+/ ^ + A . There is no such field. However, if we consider 
the density ^(aII 1 a 31 2 a 31 3 ) °^ operators when we set v4 31il = A 31j2 = A 313 = 0, then this 
density does become zero appears to be important to allow arbitrary 

sets of bifundamental fields to be set to zero. A more geometric way of saying this is that 
the important objects in the operator counting formula are not the bifundamental fields but 
rather five-cycles in the Sasaki-Einstein manifold. In the C 3 /Z 3 theory, there seems to be no 
operator constructed from bifundamental fields that corresponds to a five-brane wrapping 
the cycle A 3 i^ = A 3 i :2 = ^31,3 = 0j^] We might say that we are missing some operators. We 
note that the problem could be resolved if we added an operator A31, 1/^23,1, since the cycle 
^-31,1 = is the sum of the cycles A 23i i = and v4 311 = A 312 = v4 313 = 0. The problem 
never arises in non-chiral non-flavored theories because these theories do have an operator 
for every cycle. 

The flavored M = 8 and flavored ABJM model also have missing operators. At x = 0, 



6 Unlike setting A23.1 — ^23.2 — ^.23,3 — where there are no non- vanishing operators, when we set 
^31,1 = ^31.2 = ^31,3 = 0, the number of non-zero operators ^(^31 ± ,A 31 2 ,A 31 3 ) scales as r 3 , indicating the 
presence of a 5-cycle in the geometry. 
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there is a corner in the solutions that does not correspond to any Sy a b saturating at the R- 
charge of some bifundamental field X. Instead, the corner comes from the q fields. From the 
operator counting perspective, this corner can be explained by the fact that ip^'^ becomes 
zero at x — 0. 
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A F-Maximization for the necklace quivers 

We would like to show that to leading order in N the free energy of the necklace quivers 



with arbitrary R-charges studied in section 4.1 is maximized when i2[A tt ] = R[B a ] = 1/2 



and A m = 0. We can only show this if the gauge groups are SU(N). In the U(N) d case, 



the symmetries (2.4) imply that the free energy has flat directions, but we can nevertheless 



show that the free energy is maximized when the invariant combinations c% and C2 defined in 



eq. (4.5) are set to zero. The critical R-charges correspond to the case where there is M = 3 
supersymmetry as opposed to just Af = 2. 

The essential ingredient of the proof is the observation that the polygon V, which depends 
on c, is the polar dual of a polygon Q that does not depend on c about the unit circle centered 
at (— c/2). Let (3 a = (1, q a ) and c— (ci, C2) be vectors in M. 2 . The polygon Q is the Minkowski 
sum 

2 = {j^ Ua ^ a G R2 : Ua E ( ~ 1/2 ' 1/2) } (A-1) 
of the vectors a . Indeed, one can rewrite V as the intersection of half-planes 

V= \teM? : U-c + Yj- {uj^j < l,Wu a e (-1/2, 1/2) I . (A.2) 

The boundaries of these half-planes are precisely the polar duals of the points in Q about 
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the unit circle centered at (—c/2). 

Let Vi be the vertices of Q ordered so that the line segment between Vi and tfj+i is part 
of the boundary of Q. The line passing through Vi and Vi + \ is polar dual to a vertex w% t i+i 
of V. Polar duality implies iZ^+i ■ (ty + c/2) = w^i+i ■ {vi+i + c/2) = 1, so 

*(Vj+l ~ Vi) (A on 

M+1 " m + i + c/2)) ■ $ + c/2) ' [A - 6) 

where * denotes the Hodge dual in IR 2 . By splitting V into triangles we can write the area 
of V as 

Area("P) = ^ Area(wj_ M , w M+1 , 0) = ^ - • , (A.4) 

i i 

where we denoted the area of a triangle whose vertices are given by the vectors a, /3, and 7 



by Area(a, j3, 7). Using eq. (A. 3), eq. (A.4) becomes 



Area(-P) = - V" Area(^-i,^,^ i+ i) 

4 ^ Area(vi, -c/2) Area(v i+1 , v i} -c/2) ' ^ ' ' 

As long as —c/2 belongs to the interior of Q, the Hessian matrix of each term in this sum, seen 
as a function of c, is positive definite, so Area('P) is a convex function of c. (To compute the 
Hessian it is easiest to work in a coordinate system where c is parametrized by the distance 
from two neighboring sides of the polygon to —c/2.) 



In our case Q is symmetric about the origin as can be easily seen from eq. (A.l). Con- 
sequently, Area("P) is an even function of c, and we have just shown that it is also convex. 
It follows that Area("P) is minimized for c = 0. Equivalently, the free energy is maximized 
when c = 0. Using the F-maximization conjecture of [6], we have thus shown that the correct 



R-charges in the necklace quivers with superpotential (4.1 ) satisfy C\ = C2 = 0. That's all one 
can say about the U(N) d theory. If the gauge groups are instead SU(N), the tracelessness 
constraints J dx p(x)Sy a (x) = imply (when c = 0) 

dxp{x)8y a {x) = R[Ba] ~ R[Aa] =0, (A.6) 
so R[A a ] = R[B a ] = 1/2. From c 2 = we also get A m = 0. 
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B Further examples 

For notational convenience we set T« = T and = f. 
B.l Flavored ABJM theory 

We consider the flavored ABJM model with the superpotential 



W ~ tr 



rial 



Hal 



"61 



"62 



(2) 



i=i 



i=i 



(B.l) 



When N — 1, the superpotential is supplemented by the relation (2.7) which in this case is 



J^J^ — j^ n ai j^rial J^ribi j^ri b2 



22 



23 



The corresponding matrix model was solved in the large 
N limit in [9]. Our strategy is the same as for the flavored J\f = 8 theory In this section, 
we will review the solution for p(x) and Sy = y\ — y 2 . In the next section, we will compare 
these results with the distribution of operators in the chiral ring. 

We define R[Ai] = A Ai , R[Bj\ = A Bi , A = R[T], and A = R[f]. Without loss of 
generality, we will assume that A A2 < A Al and A B2 < A Bl . To keep the notation concise, 
we also define 



k± = k ± - (n i + n a2 - n bl - n b2 ) , 
A 2 = A Al A A2 - A Bl A B2 , 

A 3 = A Al A M (A Bl + A Ba ) + A Bl A B2 {A Al + A Aa ) . 



(B.2) 



Taking the marginality constraints on the R-charges into account, in the large N limit, 
the matrix model free energy functional is 



f\pM 

2vriV 3 / 2 



dx p 



^(k+ + k_)x 5y-p ((5y) 2 + A 2 5y- ^A 3 j + ^(A - A)x 



— pi / dx p — 1 ) . 



+ - \x\ (A + A + (k+ - k_)5y 



(B.3) 
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The eigenvalue density has four regions: 
1 1 

P : 



R[fA k 2 -} <X< R[TA\~ 



< x < 



R[fA k { 



< x < 



R[TB[ 



1 A 1 

< x < 



R[TB k 



R[TB 



1 + .Ti?[TA 2 



fc_ n 



2 ' 



A 3 + 2A A2 A 2 - 2A M 
2 + 2xA + x£;_A 2 



-Aa 2 



A| + 2A 3 

k_xA 3 - (1 + xA)A 2 



2 + 2xA + xk_A 2 
2 - 2xA + xk + A 2 
Al + 2A 3 ' 

A;+xA 3 - (1 - xA)A 2 



5£ 



2 - 2xA + xA;+A 2 ' 



A 3 - 2A B2 A 2 - 2A 



X B 2 



As in (3.2), we have introduced the rescaled variables x = x\i and p(x) = p(x)fi. 



(B.4) 



(B.5) 



(B.6) 
(B.7) 



Operator counting 

There are operators containing T~ m for m < and operators containing T m for m > 0. They 
take the form T^A^A^B^B^ 1 and f- m A?- AS? , where gauge invariance demands 
«i + a 2 — /?i — 02 = —mk±. If we wanted to count operators that don't vanish when, for 
example, A\ = 0, then we just set a± = 0. 

We counted the operators using a slightly modified version of the method outlined in 
appendix [cj Having written the operators in terms of both T and f, it is simpler to use two 
different coordinate systems on the cone C, one when m > and one when m < 0. The 



coordinate systems are related by (2.7). The operator counts reproduce (B.4), (B.5), (B.6) 



and (B.7) via our conjecture (1.4) 



Here are some of the details for the calculation of p(x) when m > 0. The density of 
operators is given by 



drdm 



da.\ da 2 d/3i df3 2 5{a\ + a 2 — (3i — (3 2 + mk^ 



x 5{r - mA - a x A Al - a 2 A A2 - (3iA Bl - /3 2 A B , 2 ) . 



(B.8) 



This integral gives the area of a slice of a tetrahedron. The slice is either a triangle or a 
quadrilateral (which may be regarded as a triangle with another triangle cut out). We find 
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h a 2 i A 12 k 2 




Figure 3: The quiver for C 3 /(Z2 x Z2). There are four U(N) gauge groups with Chern- 
Simons coefficients k a . The matter content consists of the 12 bifundamental fields A a b for 
a b, transforming under the fundamental of the 6th gauge group and the antifundamental 
of the ath gauge group. 



for m > 



drdm 



2(A 3 - 2A B .A 2 - 2AV 



(B.9) 



Taking a derivative of this expression with respect to r yields (B.6) and (B.7) 



B.2 C 3 /(^2 x Z 2 ) theory 



Let's examine the field theory in figure |3j It has four gauge groups with CS levels k a , 
a = 1, ... ,4, and twelve bifundamental fields A a b transforming in (N OJ N&), one for every 
ordered pair (a, b) with a ^ b. The superpotential is 



W = tr 



^ e abcdAdbA cc iAi JC 



a=l 



(B.IO) 



1. We 



The superpotential relations are supplemented by the monopole OPE (2.7) TT 
define R{A ab ] = R ab and R[T] = -R[f] = A. 

The superpotential contains eight distinct terms that impose the relations R a b + Rbc + 
R ca = 2 for any triplet (a, b, c) of pairwise distinct gauge groups. These eight equations 



imply the long-range force cancellation (2.1). Only seven of these equations are linearly 



independent, leaving five independent R-charges out of the twelve R a b- 

Even though for given k a the matrix model depends on 6 R-charges (A and the five 
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linearly independent R a b), the dependence on three of these parameters is trivial because of 



the flat directions (2.4). We can use these symmetries to reduce the number of independent 



R-charges to three: A x , A y and A where we pick 



Rl2 


= R21 


= R34 


= R43 


= A x , 


R23 


= R32 


= R 41 


= -R14 


= A,, 


R13 


= R31 


= -^24 


= -R42 


= 2- A x 



(B.ll) 



The matrix model is then 

d 



F[ Pl y a } = 2nN 3/2 / dx px^k a y a + 2nAN 3/2 / dx px 

J a=l J 

^^3/2 r 

H g — / dx P 2 ( yb ~U a + Rab > (V c ~Vb + Rbc) (y a -y c + R c 



(B.12) 



(a,b,c) 



For simplicity, let's focus on the case k\ = — k 2 = k 3 = —k 4 = k > and take A y > A x . 
The saddle point eigenvalue distribution splits into three regions where p is linear: 

1 - x(A - 2kA x ) 

<x< ^^- = 4A x (A y - A X )A Z ' ^-^ = - A - 

1 - xA 2kxA x A y 



1 


A 


- 2fcA x 




1 


A 


-2k A y 




1 


A 


+ 2fcA y 



1 


A 


- 2kA y 




1 


A 


+ 2k A y 




1 



<£< aT^a" : ^ ) = 4A^a? ^-^ 2 = T3I^' ( R13 ) 



1 - x(A + 2k A x ) A A 
<X< AT^ : /;! ' r) 1A,(A, A,)A/ 

In all three regions, 2/1 = y 3 and y 2 = y 4 - 



Operator counting 

Without the monopole operators, the ring of functions A a b modulo superpotential relations 
is the ring of functions on C 3 / (Z2 x Z2). This ring consists of polynomials in x, y, z with the 
constraint that the numbers of x, y, z in each term must be either all even or all odd. We call 
A 12 , A 21 , A 34 , A 43 u x fields", A u , A 41 , A 23 , A 32 u y fields", and A 13 , A 31 , A 24 , A 42 a z fields". 
We can get a gauge invariant operator by taking a combination of two x fields (e.g. A\ 2 A 2 i), 
two y fields (Ax 3 A 3 i), two z fields (Ai 4 A 4 x), or one of each type of field (Ai 2 A 23 A 3 i). The 
gauge invariant operators with m = are those with an even number of each of x, y, z, or 
an odd number of each of x, y, z. 

Adding back the monopole operators yields a ring of functions on a four-dimensional cone. 



32 




Figure 4: The area of the polygonal regions ABC and ABCD is proportional to d 2 ip/drdm 
for the C 3 /Z 2 x Z 2 quiver: a) r — Am — 2fc|m|A 2/ < 0; b) r — Am — 2/clmlA^ > 



An electric charge of (1, —1, 1, —1) from T can be cancelled out by two x's (A12A34) or two 
y's (y4 14 y4 32 ), but not by z's. So, if we have an operator of the schematic form T m x nx y ny z Uz 
for m > and T~ m x nx y Uy z nz for m < 0, then we have the constraint n x + n y > 2\m\. The 
operator density is then 



1 r 

- / rfn^ rfn^ dn z 9(n x + n y — 2k\m\)5(r — A x n x — A y n y — A z n z — Am) . (B.14) 



drdm 4 

The factor of \ comes from the constraint that the numbers of x, y, z must be all even or all 
odd. 

Performing the integral over n z introduces an overall factor of 1/A Z . The remaining 
integral reduces to the area of a polygonal region satisfying the constraints n y > 0, n x > 0, 
n x + n y > 2k\m\, and A x n x + A y n y < r — Am. For small |m|, the polygonal region is 
a quadrilateral while for large \m\, the region is a triangle (see figure [4]). Assuming that 
A y > A x , we find 



drdm 



1 (r-Am-2k\m\A x ) 2 
8A~ Z A x (A y -A x )2 



1 

8A 2 



(r— Am) 2 



(2k\m\y 



if r — Am — 2k\m\A y < , 



if r — Am — 2k\m\A y > . 



(B.15) 



Taking an additional derivative with respect to r, we can easily check that this formula 



agrees with (B.13). 



Now, in order to compute yi(x) — 2/4(2;), we count gauge invariant operators with A u 
set to zero. Because of the superpotential relations, all operators with set to zero. 

The factor of 1/4 remains the same because now we may only consider operators with even 
numbers of x and y fields. The expression for dipi^/dm is given by the area of the same 
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polygonal region that governs d 2 ip/drdm, but we lose the factor of A z because we drop the 
integral over n z : 

d 2 jj 14 = ^ d 3 jj 
drdm z dr 2 dm 

Therefore, we have p(x)(A z + yi(x) — yi(x)) = p(x)A z , and hence y\{x) = y^(x). A similar 
calculation shows ^(^) = 2/3 (^)- 

Finally we count the operators with A 12 set to zero. Most operators with an x will 
become zero. However, fields containing only T, A21, ^43, ^23, An, and the z fields are not 
set to zero by the superpotential relations. So the nonzero fields are those with n x = and 
an even number of y and z fields, or m > and n x + n y = 2km with an even number of z 
fields. After a little work, we find 



drdm 



if m/r < -(2k A y - A)- 1 

f " A g; A » if - (2kA y - A)- 1 < m/r < (2kA y + A)" 1 (B.17) 
{ r 2A^"-t) X Xrn/r>(2kA y + A)-' 



This result matches y~i(x) — y2(x) computed from (B.13). 



C Toric varieties in general 

By toric moduli space we mean more specifically that the moduli space for the Abelian gauge 
theory is an eight-dimensional toric Calabi-Yau cone V. That V is toric means it is a T 4 
torus fibration over a four-dimensional rational polyhedral cone C. This polyhedral cone is 
the set of points satisfying 

C = {yeR i :yv a >0} , (C.l) 
where v a G Z 4 , a — 1, ... ,n, are inward pointing vectors normal to the faces F a of the cone: 

F a = {y G C : y ■ v a = 0} . (C.2) 

The fact that V is Calabi-Yau implies that the end-points of the vectors v a lie in a common 
hyperplane M 3 . 

One convenient aspect of this construction is that lattice points in C correspond to 
operators in the chiral ring of the Chern-Simons theory. The coordinates of a lattice point 
are the U(l) global charges of the operator. The vector b that measures the R-charge is 
often called the Reeb vector where the R-charge is then r = y - b. The vectors v a correspond 
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to other global charges, q a = y ■ v a , and we can introduce additional charges as well. In the 
gauge theories considered in this paper, the monopole charge m played an important role. 
Let us introduce t as the vector that measures monopole charge. 

We introduced previously the function x/j(r,m) as the number of operators with R-charge 
less than r and monopole charge less than m. From the toric perspective, this function in 
the large r and m limit is the volume of a four-dimensional polytope: 

C r ,m = C n {y ■ b < r} n {y ■ t < m} , (C.3) 

where ip(r,m) = Vol(C r)TO ). 

We would like to understand geometrically how to compute derivatives of if)(r,m). The 
value of if)(r,m) is a four-dimensional integral we can write as 



ip(r, m) = 




(C.4) 



To take a derivative of ip with respect to r, we can rotate the coordinate system so that one 
of the y's points in the direction of b and replace d 4 y with d 3 y dr/\b\ where |6| is the Jacobian 
factor from the change of variables. The derivative is then related to the three-dimensional 
volume of the polyhedron 

D r , m = Cf]{yb = r}n{yt<m} (C.5) 

where dip/dr = Volp r , m )/|6|Q 

Similarly, we can visualize d 2 ip / drdm as the area of a two-dimensional polygon P r}Tn : 

P r ,m = Cn{yb = r}n{yt = m}. (C.6) 

Now we rotate our coordinate system so that two of the y's lie in the plane spanned by b 
and t. The Jacobian factor is \t A b\ = ^/t 2 b 2 — (t ■ b) 2 . Geometrically, the second partial is 

d 2 ip = Area(P r , m ) 
drdm \t A b\ 

The function ipx{ r , m ) has a toric interpretation as well. In the examples we considered, 

7 This last expression may seem strange because the right hand side seems to depend on a metric while 
the left hand side depends only on a volume form on C . Interpreting Vol(D r . m ) as a three form instead of a 
number, we could rewrite this expression in a manifestly metric independent way: (dip/dr)t = *Vol(-D rjm ). 
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X corresponds to an integer linear combination of the v a . Let us consider the simple case 
where X a corresponds to a single v a . Operators with no X a are contained in the face F a C C . 
This fact suggests a relation between ^x a (r,m) and a generalization of ifj(r,m) involving a 
third charge q a , ip(r, m, q a ). In particular, it is true that 



ip Xa (r,m) = V> (OA1) (r,ra,0) . 



Operators with no X a and fixed m and r lie along a line L a ni T C F a : 



J a,m,r 



F a n{yb = r}n{yt = m}. 



(C.9) 



Generalizing the argument used to derive (C.7) to one more charge, we find 

Length(L ai?Tlir ) 



Ml 
drdm 



V> (1,1,1) (r,m,0) 



It A 6 A U n 



(CIO) 



Eqs. (C.7) and (CIO) provide a convenient starting point for counting chiral operators in 



the examples in the text. 



D The Cone over Q 2 ^ 2 /Z k 

As another example with chiral bifundamental fields, we can examine the square quiver in 
figure [5] with CS levels (k,k,—k,—k) and matter fields A4, Bi, Ci, and D^, with i — 1,2. 
With the superpotential is 



W ~ tr [e^e kl D i C k B j A l ] 



(D.l) 



this quiver is thought to be dual to AdS^ x Q 2 ' 2,2 /Z* 29 30 . The quiver has two flavor 
SU(2) symmetries, one under which Ai and Cj transform as doublets, and one under which 
Bi and transform as doublets, so one expects the R-charges of the fields belonging to the 



same edge of the quiver to be equal when F is maximized. Using the flat directions (2.4) 



and taking into account the marginality of the superpotential (D.l), one can then set the 
R-charges of all the bifundamental fields equal to 1/2 and A m = 0. With this choice one can 
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A; 



Di 



-k 



C; 



B: 



-k 



Figure 5: Quiver gauge theory believed to be dual to AdS± x Q 2 ' 2 ' 2 /^ 



50 through the operator counting exercise in the Abelian theory and predict that 

p{x) = 



1 

2k 



1 -1 l 1 r ( 1 A 1 r f 1 

\x\ H — -0 — + a; H — -0 — - — x 



(y 2 (x) - yi(x)) = --S 



8k 6 \2k X 



p{x) (ys{x) - y 2 {x)) 



As a consistency check, one can compute the volumes 



Vol(F) 

Vo1(SaJ 
Vol(E Bi ) 
Vol(S c J 
Vol(S A ) 



71 

24 

T 

4 
4 

4 



(if 



7T 



16fc 

rfxp(x) ( y 2 (a) - yx(x) 



dxp(x) \ y 3 (x) - y 2 (x) + 
dxp(x) H)4(x) - + 
dxp{x) ( yi(£) - £4 (£) + 



7T 

4fc ' 
8A; ' 
4A; ' 



(D.2) 



(D.3) 



Since Vol(Q 2 ' 2 - 2 ) = vr 4 /16 25 , we see that Vol(F) matches that of a Z k orbifold of Q 2 > 2 ' 2 . As 



for M 1,1 ' 1 , we can relate the volumes of the five-cycles in (D.3) to those computed in [25]. 



The cone over Q 2,2 ' 2 is a U(l) 2 Kahler quotient of C 6 with weights (1,1,-1,-1,0,0) and 
(1,1,0,0,-1,-1), together with a Z 2 quotient that flips the sign of (ai,a 2 ). If we denote 
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the coordinates in C 6 by (ai, a 2 , £>i, &2, Ci, c 2 ), we have [25] 

Vol(g 2 ' 2 ' 2 ) = ^ , Vol(E a J = Vol(EO = Vol(E c J = £ . (D.4) 

One can think of the Z^ orbifold as acting on Cj with opposite phases, so it is natural to 
interpret Ai and Cj as correponding to a i; i?j as corresponding to fejCi, and Dj as correspond- 
ing to b iC2 . Indeed fcVol(E A ) = Vol(E c J = Vol(E a .), fcVol(E fli ) = Vol(E 6i ) + Vol(E cl ), and 



k Vo^E^J = Vol(E6j + Vol(E C2 ), the factor of k appearing because the volumes (D.3) are 
computed in a Z& orbifold of Q 2,2 ' 2 . 
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